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i . For the one-dimensional XXX model under the periodic boundary condi- 

tions, we discuss two types of eigenvectors, regular eigenvectors which have 

o ; 

finite- valued rapidities satisfying the Bethe ansatz equations, and non-regular 

cn . 

J> . eigenvectors which are descendants of some regular eigenvectors under the 

o 

^O : action of the SU(2) spin-lowering operator. It was pointed out by many 

authors that the non-regular eigenvectors should correspond to the Bethe 



ansatz wavefunctions which have multiple infinite rapidities. However, it has 
not been explicitly shown whether such a delicate limiting procedure should 
be possible. In this paper, we discuss it explicitly in the level of wavefunc- 



a. 

tions: we prove that any non-regular eigenvector of the XXX model is derived 
O 

from the Bethe ansatz wavefunctions through some limit of infinite rapidities. 
We formulate the regularization also in terms of the algebraic Bethe ansatz 
method. As an application of infinite rapidity, we discuss the period of the 
spectral flow under the twisted periodic boundary conditions. 
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I. INTRODUCTION 

The one-dimensional Heisenberg model (XXX model) under the periodic boundary con- 
ditions is one of the fundamental models of integrable quantum spin systems. |lj Under 
the spin SU(2) symmetry any eigenvector of the Hamiltonian is given by a highest-weight 
vector or a descendant of some highest-weight vector. It has been shown by the algebraic 
Bethe ansatz method that any regular Bethe ansatz eigenstate of the XXX model gives 
a highest weight vector. 0,[§]. Let us consider the XXX Hamiltonian under the periodic 

boundary conditions 

J L 

TC = -— ^2&e ■ &i+i, where a L+1 = B x . (1.1) 

4 i=\ 
Here the symbol oi = (af , a\ ', of) denotes the spin angular-momentum operator with S = 

1/2 acting on the £-th site of the ring. Let us denote by the symbol S to t the total spin- 
operator: Stot = Y,£=i &e/2- Then, it is easy to show that the Hamiltonian is invariant under 
the action of the SU(2): [H, S tot ] = 0. 

Let us introduce some notation of the coordinate Bethe ansatz [[|,[5]J|. We denote by 
x\,X2, ■ ■ ., xm the coordinates of the M down-spins set in increasing order: 1 < X\ < 22 < 
• • • < Xm < L. Then, we define the Bethe ansatz wavefunction with M parameters ki, k 2 , 
. . ., ku by the following: 

/ M \ 

/If (^l,- ..,x M ;ki,...,k M ) = Yl A M (P)expliJ2kpjXj\, (1.2) 

PeS M V J =1 / 

where the sum is over all the permutations of M letters of the set {1,2,..., M}, and the 
symbol Pj denotes the action of permutation P on letter j . Here the symbol Sm denotes 
the permutation group of M letters. We define the amplitudes Am (-P)'s of the Bethe ansatz 
wavefunction by 

A M (P) = Ce(P) n ^l'(^ + ^)] + 1 -^P^) , for PeSM . (1 . 3) 

i<j<e< M ex PK fc i + ki)\ + 1-2 exp(ikj) 

Here the symbol e(P) denotes the sign of permutation P, and C is a constant. Let the 
symbol |0) denote the vacuum state where all spins are up (M = 0). Then, we construct 
the following vector from the Bethe ansatz wavefunction 



\\M)= £ fif\x lj ... J x M ;k l ,...,k M )a- i a- 2 ...a- M \0). (1.4) 

l<xi<a:2<*"<a;jwr <L 

Here, the summation is over all the possible values of x/s given in increasing order. We call 
the vector \\M) ( |1.4|) with the amplitudes defined by eq. (|1.3|) , a formal Bethe vector (or 
formal Bethe state). We recall that there is no constraint on the M parameters k\, fe, • • •, 
ku- When they are generic, the formal Bethe state (|1.4| ) is not an eigenvector of the XXX 
Hamiltonian. 

Now, let us consider the Bethe ansatz equations. They correspond to the periodic bound- 
ary conditions for the Bethe ansatz wavefunction. 

cti \ i 1 NAf-i tt ex P[#j + h)] + 1-2 exp(ifcj-) 

exr>[iLki) = ( — 1) ttt^ r4i t-tt, 

yK 3) K ' e= \ i i _ tj exv[i{k j + k i )] + l-2eMih) 

for j = l,...,M. (1.5) 

If all the parameters ki, ki, ■ ■ ., ku satisfy the Bethe ansatz equations, then the formal Bethe 
vector \\M) becomes an eigenvector of the XXX Hamiltonian. Furthermore, if the k/s satisfy 
the conditions that kj ^ (mod27r) for j = 1, . . . , M, then we call the eigenvector regular, 
and denote it by the symbol \M). It is called regular, since it is well defined as an eigenstate 
given by the Bethe ansatz wavefunction. In this sense, it is also called a regular Bethe ansatz 
state or a Bethe state, in short. 

A regular eigenstate can lead to a series of non-highest weight eigenvectors of the SU(2) 
symmetry. Let \R) denote a given regular eigenstate with R down-spins. Then, it is a 
highest weight vector of the SU(2) symmetry with S tot = L/2 — R and S* ot = L/2 — R. 
Here we assume that the number R should satisfy the condition: < R < L/2, for regular 
eigenvectors. From the eigenvector \R), we can derive a sequence of non-highest weight 
eigenvectors: (Sj~ ot ) K \R) for K = 1, . . . , L — 2R. We call the series of descendant eigenstates 
non-regular. We denote them by 

\R,K) = ^(S M ) K \R) for K = 1,...,L-2R. (1.6) 

It is remarked that the eigenvectors \R, if)'s are fundamental in the completeness of the 
spectrum of the XXX model, although they are called non-regular in this paper. 



The main question of this paper is how non-regular eigenvectors of the XXX model are 
related to the Bethe ansatz wavefunctions. In fact, it has already been observed by Gaudin 
jTj that the non-regular eigenvectors are associated with the Bethe ansatz wavefunction 
with several parameters k/s being equal to zero. Furthermore, it was shown by Takhtajan 
and Faddeev || that the creation operator B(v) is equivalent to the spin- lowering operator 
Stat by sending the rapidity v to infinity (see also Refs. [[SHTTll). We note that the limit of 
sending the parameter kj to zero corresponds to the limit of the infinite rapidity. Here, for a 
given parameter k, the rapidity v has been defined by the relation: exp(ik) = (v + i)/(v — i)] 
rapidity v is finite if and only if k ^ (mod 2tc) . In spite of the observations, however, it has 
not been clearly shown yet whether one can construct the non-regular eigenvector \R, K) 
from the Bethe ansatz wavefunctions for the case of general K. In the case of multiple 
infinite rapidities, the limit of the wavefunction depends not only on its normalization but 
also on how we control the differences among the infinite rapidities. Thus, under a naive 
limiting procedure, the amplitudes of the formal Bethe state become indefinite; it can vanish 
or diverge depending on the limiting procedure. (For example, see also \T2\.) Furthermore, if 



a set of parameters kj's contains multiple zeros, then it is not clear whether the Bethe ansatz 
wavefunction should vanish or not. In fact, for any given regular eigenvector, we can show 
that if two momenta (or two rapidities) have the same value, then the norm of the eigenvector 
is given by zero. This fact is called the "Pauli principle" of the Bethe ansatz wavefunction. 
Thus, the question has been nontrivial. In this paper, we make it clear. We show that there 
exists a certain limiting procedure through which any non-regular eigenvector of the XXX 
model is derived from the formal Bethe state. 

Let us explain our derivation of non-regular eigenvectors from the formal Bethe states, 
briefly. We consider a given regular Bethe ansatz eigenstate \R) with R down-spins. It has 
R rapidities Vi,v 2 , ■ ■ ■ ,Vr, satisfying the Bethe ansatz equations for R down-spins. For a 
given positive integer K, we consider the non-regular eigenstate \R, K). We recall that it 
has been defined in eq. fll.GQ and is derived from \R). Then, we introduce an additional set 
of the rapidities Vr+i, ■ ■ ■, vr+k as follows 



v R+j (A)=A + 6 j , for j = l,...,K. (1.7) 

Here we call the parameter A the "center" of the additional K rapidities Vr+i, ■ ■ ■ ,vr+k- 
We assume that the 5j 's are arbitrary non-zero parameters, which can be sent to infinity. 
Let us now consider a formal Bethe vector \\R + K) with R + K down-spins that has R 
rapidities of the given regular eigenstate \R) (i.e., vi, . . . ,vr) together with the additional 
K rapidities given by eq. ( |1.7|) (i.e., vr + i(A), . . . ,vr + k(A)). We denote it by \\R,K;A). 
Then, we can show that the vector \\R,K; A) becomes the non-regular eigenstate \R, K) by 
sending A to infinity: 

lim \\R,K;A) =C\R,K) (1.8) 

A^oo 

Here C denotes a constant. Thus, the non- regular eigenstate is derived from the Bethe 
ansatz wavefunction. 

We discuss only regular eigenvectors of the XXX model and their descendants which we 
call non-regular eigenvectors. We do not consider other types of solutions in this paper. 
In fact, it was shown that the so-called string hypothesis predicts the correct number of 
appropriate solutions to the Bethe ansatz equations of the XXX model under the periodic 
boundary conditions lHJD^pfl . Although the hypothesis fails to count the particular type 



of solutions, all the known numerical or analytical researches have shown that the total 
number of solutions to the Bethe ansatz equations is given correctly |I| Jl5| 1 |l6| . |TT| . Thus, it 



is conjectured that all the regular eigenvectors and their descendants give the complete set 
of eigenvectors of the XXX model. In fact, it is proven that the number of solutions of the 
Bethe ansatz equations is given correctly for the XXX model under the twisted boundary 



conditions with the generic twisting parameter [Tj. It seems that the theorem does not 
cover the case of the periodic boundary conditions, since it corresponds to a non-generic 
point of the twisting parameter. However, the result of the paper might also shed some 
light on the mathematical understanding of the string hypothesis and the number counting 
arguments in general, as we shall discuss in sections 5 and 6. 



The contents of the paper consists of the following. In section 2 we give a formula 
describing the action of powers of the spin-lowering operator. Then, through some examples, 
we explicitly discuss the derivation of non-regular eigenvectors from formal Bethe states. It 
is shown that infinite rapidities do not always satisfy the Bethe ansatz equations, although 
the limit of the Bethe ansatz wavefunction satisfies the periodic boundary conditions. In 
section 3, we give an explicit proof for the construction of non-regular eigenstates from 
the formal Bethe states. We show in section 4 that the formal Bethe state can be defined 
naturally in the algebraic Bethe ansatz. In fact, the formal Bethe state \\M) is equivalent 
to the vector generated by the B operators on the vacuum: B(v{) ■ ■ ■ B(vm)\Q) with the M 
rapidities V\, . . . ,vm being generic. In section 5, we show how the concept of formal Bethe 
states is useful in the analysis of the spectral flow of the XXX model under the twisted 
boundary conditions. In fact, we can derive the 47r-period of the spectral flow under the 
twisted boundary conditions, almost rigorously. In section 6, we give some discussions. In 
order to make the paper self-consistent, some Appendices are provided. The formula for 
the action of spin-lowering operator is proven in Appendix A. An example of the formal 
Bethe state with three infinite rapidities is discussed in Appendix B. Some fundamental 
properties of the symmetric group are given in Appendix C, which are important in sec. 3. 
The "Pauli principle" of the Bethe ansatz wavefunction is explicitly proven in Appendix D. 
Finally, we formulate rigorously the coordinate Bethe ansatz method introduced by Bethe 
|J in Appendix E. 

II. FORMAL BETHE STATES AND NON-REGULAR EIGENSTATES 

A. Non-regular eigenstates 

Let us discuss the action of spin-lowering operator on arbitrary vectors with M down- 
spins, explicitly. For an illustration we consider the case of M = 1. Let |1) denote a vector 
with one down-spin 



|i)= Z^iK~l°>> (2- 1 ) 

£1 = 1 



where g(x) is any given arbitrary function. Applying to it the spin- lowering operator S, 



tot 



Ey=i a j j we nave 



^=1 3 

L L 

s tot\ l ) = Z ^ Z ^O^J ) 

12=1 XI 

L L 
= Z Z ^OO^i^l ) 

£1 = 1 £2 = 1 

Z + Z WiK^Jo) 

,l<a;i<a;2<i l<x2<xi<Lj 

Z (^O + ^O^K^Jo) 

l<xi<£2<£ 

= Z (Z^i))^^l°) (2-2) 

l<X!<x 2 <L \1<J<2 / 

Here we note that (<r~) |0) = 0. 

We can generalize the expression (|2.2|). Let us denote by the symbol \M) a vector with 
M down-spins 

\M)= J2 9(x 1 ,x 2 ,---,x M )(r n (r^---(r- M \0), (2.3) 

1<xi<X2<---<xm<L 

where g(x±,X2, ■ • • ,xm) is an arbitrary function of x/s. Then, it is clear that any vector 
with M down-spins can be considered as a vector \M) with some function g(xi, x 2 , ■ ■ • , %). 
Now, we introduce the following formula 



-^{S M ) K \M)= Z ( Z g(x n ,... } x ni ))a^...a~ M+K \0). (2A) 

1<X1<-<X M +K<L \l<jl<-<JM<M+K J 

We note that the expression (|2.2|) corresponds to the case M = K = 1. An explicit proof 



of the formula (2.4) will be given in Appendix A. In sec. 2.C, we shall consider the special 
case of K = 2 and M — 1, which is given in the following 

1 2 

2^) I 1 ) = Z Z dix^a'^a-jO) 

"~ l<xi<X2<a;3<L 1<J<3 

Z {9(x 1 )+g(x 2 )+g(x 3 ))(T- 1 (T^a- 3 \Q). (2.5) 

l<xi<X2<a;3<L 



Here we note that M + K = 1 + 2 = 3. 

Let us consider a regular eigenstate \R) with R down-spins, and the non- regular eigen- 
state \R, K) given by eq. ( |1.6| ). We recall that \R) is a highest weight vector of the 577(2) 
with 5* = L/2 — R and S z = L/2 — R. By applying the formula fl2.4| ) to the definition ( |1.6| ) 
of the non-regular eigenvector, then it is explicitly expressed in terms of the Bethe ansatz 
wavef unctions 

\R,K)= E ( E /i?W--.^W"Wl°>- (2.6) 

l<x 1 <-<x R+K <L \l<ji<-<JR<R+K j 

Here we recall that the function f R (xi, ■ ■ ■ xr] hi, . . . , Um) is the Bethe ansatz wavefunction 
defined in eq. (|1.2| ), where the kj 's satisfy the Bethe ansatz equations. 

B. Amplitudes of formal Bethe states 

Let us recall the relation between rapidity Vj and parameter kj 

exp{ikj) = ^A for j = 1, . . . , M. (2.7) 

In terms of rapidities, the Bethe ansatz equations are given by 

(*±4V= fl f^±f), for i = l,...,M. (2.8) 



The amplitudes Ajy(P)'s defined in eq. ( p..3|) are given by the following 



A M (P)[v u ...,v M }=e(P) n 7_r+ + 2 f (2 - 9) 

Here, the dependence of the amplitude Ajvf(P) on rapidities fi, . . . ,% is expressed in the 
bracket [••■], explicitly. Here we note that the expression ( |L~3| ) of the amplitude A M (P) is 
proven in Appendix E. 

Let us now introduce a useful formula for expressing the amplitudes of the Bethe ansatz 
wavefunction. We denote by the symbol H(x) the Heaviside step function defined by H(x) = 
1 for x > 0, and H(x) = otherwise. Then, we can show that the amplitudes Am(P)'s 
given in eq. (|2.9| ) are expressed by the following 

8 



Vi — v 2 + 2iJ \vi - v 3 + 2i 

Vi — v 3 — 2i\ fv 2 — v 3 — 2i 



/ -\ H(p- 1 j-P~ 1 k) 

l<i<fc<A/ Vo ^ + Z V 

We shall prove the expression (|2.10|) in sec. 3. 

For an illustration, we consider the amplitudes Am(P)'s for the case M — 3. Let us 
express Am(P) by Ap\p 2 ...pM- Then, they are given as follows 

a _i a _v 2 -v 3 -2i _v 1 -v 2 -2i 

^123 — -L; ^132 — — pTT, ^213 — ~ TT, 

^2 - ^3 + 2l Vl — V 2 + 2l 

(V\ — V 2 — 2i\ (V\ — V 3 — 2i' 
^231 

312 

Vi — Va + 2iJ \v 2 - v 3 + 2% 

A 321 = f v ^-^- 2 ') (*-*-*) (vi-v*-*) (2 . n) 

\Vi — v 2 + 2iJ \vi — v 3 + 2i/ \v 2 — v 3 + 2z 



C. Formal Bethe states with additional infinite rapidities 

Let us discuss some examples of the Bethe ansatz wavefunctions with additional ra- 
pidities. We first consider the case of three down-spins with R = 1 and K = 2, i.e., the 
formal Bethe state ||1,2; A). Here, v 2 and v 3 are additional rapidities defined by eq. ( |1.7D : 
v 2 = A + 6i, v 3 = A + 5 2 . Here we assume that Si and 5 2 are some constants. We recall that 
Vi is the rapidity of the state |1) and it satisfies the Bethe ansatz equation for M — 1. 

Let us denote the difference 5\ — 5 2 by A. For simplicity, we assume that 5\ = —S 2 . 
Then, the additional rapidities are given by v 2 = A + A/2 and v 3 = A — A/2. Substituting 
the rapidities vi, v 2 and v 3 into the amplitudes in ( j2.11| ), we have 



a /an 1 a /an A - 2z , , . . Vl - A - A/2 - 2i 

A 123 (A) = 1, A 132 (A) = , A 2 i 3 (A) = — ; -!-, . 

123V ; ' 132V ; A + 2*' " v ; Vi-A- A/2 + 2i- 

'v x - A - A/2 - 2i\ fv 1 -A+ A/2 - 2z\ 



^231 (A) 
^312 (A) 



y x - A - A/2 + 2%) \vi - A + A/2 + 2i / 
'vi - A + A/2 - 2z\ /A - 2z s 



vi - A + A/2 + 2%) \A + 2z 



/ % - A - A/2 - 2A U - A + A/2 - 2i \ (A - 2i\ 



Let us denote by f^ K the Bethe ansatz wavefunction for the formal state \\R,K;A). The 
Bethe ansatz wavefunction of ||1, 2; A) is given by 

fi B 2 \x 1 ,x 2j x 3 ;k 1} fc 2 (A), k 3 (A)) 
= A 123 exp i (hxi + k 2 (A)x 2 + k 3 (A)x 3 ) + A l32 exp i (k i x l + k 3 (A)x 2 + k 2 (A)x 3 ) 
+ A 213 exp i (k 2 (A)x l + k x x 2 + k 3 (A)x 3 ) + A 312 exp i (k 3 (A)xi + k x x 2 + k 2 (A)x 3 ) 
+ A 231 exp i (k 2 (A)xi + k 3 (A)x 2 + hx 3 ) + A 321 exp i (k 3 (A)xi + k 2 (A)x 2 + kix 3 ) 

for 1 < xi < x 2 < x 3 < L , (2.13) 

where k 2 (A) and k 3 (A) are given by 

-™=(£^)> -«A))^(^|±|). (2,4) 

Sending the center A to infinity : A — > oo, we have &2 = &3 = (mod)27r and 

^123(00) = ^213(00) = ^231(00) = 1 

A — 2i 
A 132 {cx>) = A 312 (oo) = A 321 {oo) = - . (2.15) 

A + 1% 

Therefore, the limit of the Bethe ansatz wavefunction is given by 

lim fi B 2 \x 1 ,x 2 ,x 3 ;k 1 ,k 2 (A),k 3 (A))=C 2 ( e <*i*i + e *i»* + e *i*s) . (2.16) 

where the constant C2 is given by 



Combining the eqs. (|2.16|) and (|2.5| ), we obtain the following result. 



A lim ||1, 2; A) = C 2 £ (e** 1 + e lk ^ + e ik ™) er" <r" <r" |0) 

l<Xi<X2<X3<L 
= ^2^(^) 2 |1) 

= C 2 |l, 2) (2.18) 

Thus, we have shown that the limit of the formal Bethe state ||1,2; A) is equivalent to the 
non-regular eigenstate |1,2). We shall prove this equivalence for the general case in sec. 3. 
For an illustration, we shall consider the case of R = and K = 3 in Appendix B. 

10 



Let us give some remarks on eq. ( |2.18|) . We see that the limiting procedure depends on 



the difference A. If A = — 2i, then the constant C<i becomes infinite. If A = 0, then the 
constant C^ vanishes. Thus, the limit of the wavefunction with infinite rapidities v 2 and v 3 
depends on how we send them into infinity. 

D. The P.B.C.s for the limits of the formal Bethe states 

The formal Bethe state | \R, K; A) satisfies the periodic boundary conditions after taking 
the limit: A — > 00. In fact, it is clear since the limit gives the non-regular eigenvector 
\R,K), which satisfies the periodic boundary conditions. Here we note that the total spin 
operator Stu is translation invariant. However, infinite rapidities do not always satisfy the 
Bethe ansatz equations. 

For an illustration, let us consider the formal Bethe state ||1,2;A). We denote by 
/i,2 {x\,X2,x^) the limit of / 12 \xi,X2,x 3 ] k\, & 2 (A), ^3 (A)) with A sent to infinity. We see 
that it satisfies the periodic boundary conditions: f\^ (x±, X2, X3) = fi°2 ( x 2,Xz,Xi + L) for 
1 < Xi < x 2 < x 3 < L. Explicitly we have 

9 A 

fif(x 2 ,x 3 ,x 1 + L) = ^-p^ (e ifclX2 + e ik ^ + e ik ^ +L ^ . (2.19) 

Thus, it satisfies the periodic boundary conditions if and only if the following holds: 

exp(iAaL) = 1 (2.20) 

This is nothing but the Bethe ansatz equation for fci, and it does hold from the assumption 
that V\ is the rapidity of a regular eigenvector |1). 

Now let us show that the additional rapidities do not necessarily satisfy the Bethe ansatz 
equations, although the limiting Bethe ansatz wavefunction satisfies the periodic boundary 
conditions. Let us consider the Bethe ansatz equations for three rapidities i>i,i>2 and V3 

v 1 + i \ L ( v 1 — v 2 + 2i \ ( v\ — v 3 + 2% 



V\ — i) \vi — V2 — 2i J \vi — v 3 — 2i 

v 2 + i\ L fv 2 -vi + 2i\fv2~v 3 + 2i 



V2 — i/ \V2 — v\ — 2iJ \V2 — V3 — 2i 

11 



V3 + i \ L = / v 3 -vi + 2i \ / v 3 -v 2 + 2i \ 
sV 3 -i) \v 3 —vi - 2%) \v 3 - v 2 - 2%) 

Taking the limit: A — ► oo, the three equations are reduced into the following 

-^— = 1 2.22 

Vx — lJ 

The equation ( 2.23| ) does not hold if A takes a finite value; it holds only if |A| = oo. 



III. PROOF OF THE LIMIT OF FORMAL BETHE STATES 

In this section we prove the theorem in the following. 

Theorem III.l Let \R) be a regular Bethe ansatz eigenstate with R down-spins and ra- 
pidities v i, . . . , vr. We recall that the symbol \\R,K; A) denotes the formal Bethe state with 
R + K down-spins, which has the R rapidities v\, . . . , Vr of \R) together with additional 
rapidities vr + x(A) . . . ,vr + k(A)- Then, the non-regular eigenstate \R,K), which is a de- 
scendant of R, is equivalent to the limit of the formal Bethe state \\R,K;A) with A sent to 
infinity: 

lim \\R,K;A} = C K \R,K) (3.1) 

A^oo 

A. Derivation of the formula for amplitudes Ajvf(P)'s 

We now discuss the derivation of the formula (|2.10|) , which rewrites the amplitudes 
Am(-P)'s defined in (|1.3|). Let us recall that the symbol H(x) denote the Heaviside step 
function defined by H(x) = 1 for x > 0, and H(x) = otherwise. We now show the 
following proposition. 

Lemma III.l Let P be an element of Sm, and v i, t>2, . . . , vm be generic parameters. Then, 
the following identity holds. 

12 



(3.3) 



(\ H(P~ 1 j~P~ 1 k) 
v k - Vj + 2i \ 
Vj-v k + 2i) 

(Proof) Let us take a pair of integers j and k with j < k, and consider the factor Vj — Vk + 2i 
in the denominator of LHS of eq. (|3.2| ). For the pair, there exist two integers £ and m such 
that P£ = j, Pm = k. There are two cases either £ < m or I > m. If £ < m, then we 
have the factor vpi — vp m + 2z in the enumerator of LHS of eq. ( |3.2|) . Thus, the factors 
associated with the rapidities Vj and v k cancel each other. On the other hand, if £ > m, we 
have vp m — Vp# + 2i in the enumerator of LHS of eq. ( |3.2| ), and we have 

VPm — Vp£ + 2i Vk — Vj + 2% 

Vj — Vk + 2i Vj — Vk + 2i 
We can express these results by the following 

/ , -\ H(i-m) 

( Vk- Vj + 2t \ 
\Vj -v k + 2i) 

Considering all the pairs j, k with j < k, we establish the equality (^ 

Q.E.D. 



Proposition III.l The amplitude Am{P) defined by eq. ( \1 . 3j ) for P G Sm can be expressed 
as follows. 

/ n .\ H(P- 1 j-P- 1 k) 

A ^- n fe^f (3.4) 

l<j<k<M \ U 3 Uk ^ Z V 

(Proof) The amplitude A M (P) defined by eq. ( |1.3|) is written in terms of rapidities as follows 

l<j<k<M U 3 Uk "^ ^ L 

In Appendix C, we show the following identity in Prop. CI 

cm{P)= II (-l) H ( p - lj - p - 1 V (3.6) 

l<j<k<M 

Thus, making use of Lem. III.l and Prop. CI, we obtain 

13 



H(p- 1 j-P- 1 k) 

— Vj T Z,l I 
l<j<k<M \ u 3 ~ 



l<j<k<M V J Vk + ^ l 

^ n terfY 



- n fefrl ■ (3.7) 

We give a remark. Using Prop. III. 1, we can explicitly prove that the Bethe states (and 
also the formal Bethe states) should vanish when there are two momenta of the same value. 
The proof is given in Appendix D. 

B. Proof of the limit 

Let us take a permutation P on R + K letters (P G S r+ k)- We consider the following 

set 

P- 1 {1,2,...,R} = {P- I j | for j = l,2,...,R}. (3.8) 

Let us denote the elements of the set by ai, a 2 , ■ ■ ■ , a R , where Oj's are set in increasing order: 
a\ < a 2 < ■ ■ ■ < a R . For the permutation P, we introduce permutation P R on R letters by 

P R m = Pa m for m — 1, . . . , R. (3.9) 

Then, we have the following. 

Lemma III. 2 Let P R denote the permutation on R letters defined by ($■§ ) for a given 
permutation P on R + K letters. For two integers j\ and j 2 with 1 < jx,j 2 < R> the 
inequality P~ 1 ji < P~ l j2 holds if and only if Pr~ ji < Pr~ J2- Equivalently, we have 

H(p- 1 3l -P- 1 j 2 ) = H(P R - 1 j 1 -P R - 1 j 2 ) for l<jt,j 2 <R. (3.10) 
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(proof) Let us denote P~ 1 ji and P~ 1 j 2 by a mi and a m2 , respectively. Then, by definition, 
we have m-i = P R ~ 1 ji and m 2 = Pr~ J2- Here we recall that a/s are set in increasing order. 
Thus, we see that o mi < a m2 if and only if mi < 1712, which gives the proof. 

Similarly, let us introduce a permutation on K letters. We consider the following set 

P- 1 {R+1,R + 2,...,R + K} = {P- 1 J | for j = R + 1, R + 2, . . . , R + K}. (3.11) 

We denote by 61,62, • • • ,bx, the elements of the above set. Here we assume that b/s are 
in increasing order: 61 < 62 < • • ■ < bx- We define permutation Pk on K letters by the 
following 

P K m = Pb m for m = l,2,...,K. (3.12) 

Then, we can show the following. 



Lemma III. 3 Let Pk denote the permutation on K letters defined by ( 3.13j ) for a given 



permutation P on R + K letters. For two integers j± and j 2 with R + l < ji,j 2 < R + K, the 

inequality P~ x ji < P~ x ji holds if and only if Pk~ (ji — R) < Pr 1 (J2 — K). Equivalently, 

we have 
H(p- 1 Jl -P- 1 j 2 ) = H(P K - 1 (j 1 -R)-P K - 1 (j 2 -R)) for R+ 1 < j 1: j 2 < R+ K. 

(3.13) 

Making use of Lemmas III.l, III. 2 and III. 3, we now show the following proposition 

Proposition III. 2 Let us consider two positive integers R and K satisfying < K < L — 
2R. Let Vi, v 2 , ■ ■ ■ , Vr be the rapidities of a given regular eigenvector \R) with R down-spins, 
and Vr + i(A), . . . ,t> fi+ ^(A) be additional K rapidities which are given by v R+ j(A) = A + 5j 
for j = 1,2, ... ,K. Here Sj 's are arbitrary constants. For the Bethe ansatz wavefunction 
fn+K with its amplitudes A R+K (P) ? s given by (\l.ty , we have the following limit. 

lim f R+K (xi, ..., x R+K ; fa, . . . , k R , k R+1 (A), . . . , k R+K (A)) 
= C K Yl fn( x jn ■ ■ ■ > x jr,i k u---i k R ) (3.14) 

l<h<-<3R<R+K 
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Here kj 's are related to the rapidities Vj 's through the relation: expikj = (vj + i)/(vj — i), 
and the constant Ck is given by 

C K = J2 A K (P)[6i,---,6 K ] (3.15) 

Pes K 

(Proof) We recall that the Bethe ansatz wavefunction f R+ x is given by 

R+K 

f(x 1 ,...,x R+K ) = J2 A R+K (P)exp(i J2 k PjXj) (3-16) 

Pes R+K j=i 

Let us take a permutation P in S r+ k- By Lemma III. 2 we can show that the amplitude 
Ar+k{P) of the formal Bethe state is given by 

/ _ _ 0'\ H(P~ 1 j-P- 1 i) 

A R+K (P)[v 1 ,...,v R ,v R+1 (A),...,v R+K (V} = II \ Vi ~ Ve lo % ) (3 - 17) 

l<j<i<R+K \Vj-Vi + 2tJ 



The above product can be decomposed into the three parts in the following 

/ Vj - Ve _2i \ = 17 ( v i- v *- 2l \ 

i<j<i<R+K \ v i -ve + 2i) i< 3 <t<R \ v i -v t + 2i) 



/ n .\ H(P- 1 j-P~ 1 i) 

x n n F^ri 

l<j<R R+l<l<R+K \ u i Vi^^l) 

/ r>-\ H(P- 1 J-P~ 1 1) 

x n 



R+l<j<i<R+K \ V J Vi + 2% t 



(3.18) 



First, we consider the third part of RHS of ( 3. 18] ). Making use of Lemma III. 3, we have 



/ n.\H(P- X 3-P-H) / n.sHiP^j-P-H) 

yr f V j -V e -2l \ = 17 \ VJ+R ~ V£+R I 

R+l<j<J<R+K \ V J -V t + 2i) 1<3<\<K \ V J+R ~ V t+R + 2 V 

=iAAi^T2i) (3 - 19) 

We note that RHS of ( |3.19| ) is nothing but A K (P K ) [Si,- ■ ■ ,5k}- Second, it is clear that 



the second part of RHS of (|3.18 ) becomes 1 under the limit: A — > oo. In fact, putting the 



additional rapidities into the second part of RHS of ( p,18|) , we have 



/ n .\ H(P- 1 j-P- 1 t) 

n 



i<j<R r+i<£<r+k \ v i ve + 2t / 

n n ( Vi z A *~l e '7- ) ■ ( 3 - 2 °) 

'-j<R R+KKR+K \ V J A de + ZlJ 



a e- «-\ H(P- 1 j-P- x tj 

Vj- k-8i-2i\ 



l<j<R R+1<£<R+K 
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Third, we consider the first part of RHS of ( |3.18|) . We recall that Pr is defined for the given 



permutation P by the relation ( |3.9| ). Then, from Lemma III. 2, we have 

n h- v <-*) = n h- v <-*) (3.2i) 

l<3<i<R \ V i ~ V * + 2 V l<j<i<R \ V J -V£ + 2iJ 



We note again that RHS of ( |3.21| ) is equal to A R (P R ) [vi, • ■ ■ , vr]. Thus, we have 



lim A P+K {P) [v 1 , ■ ■ ■ , v R , v r+1 (A), ■■■, v R+K (A)} 

A^oo 

A R (P R )[v ir --,v R ] x A K (P K )[6 U ---,8 K ] (3.22) 



Let us now consider the exponential part of (|3.16|) . We note the following 



R+K R+K 

Y^ kpjXj = Y hxp-H 

R R+K 

= J2k e x P ~i e + J2 hxp-u- (3.23) 

l=\ l=R+\ 

Since k R+ i . . . k R+ x are approaching to (mod 2tt) in the limit of sending A to infinity, we 
have 

R+K 

lim V k e (A)x P -it = (mod 2tt). (3.24) 

A ^°°e=R+i 

Making use of the relation P R m = Pa m , we have 

R R R 

2 , kiXp-i£ = 2^ kp Rm xp-ip Rm = 2_^ kp Rm Xa m (3.25) 

1=1 m=l m=l 

Thus, we have 

R+K 

lim A R+K (P)exp( V k Pj Xj) 

A — >oo . * 

R 

= A K {P K )[5 1 ,...,5 K }xA R (PR)[v 1 ,...,v R \exp{Y,kp Rm x am ) for P e S R+K (3.26) 

m=l 

Finally, we give a remark. To pick up a permutation P on R + K letters is equivalent to 
do the procedures in the following: we take a subset {ax, 02, • • • , a R } of the set of i? + -K" 
letters 1, 2, . . ., R + K, and specify Pr on i? letters and Pr on if letters by ( |3"D| ) and ( |3.12| ), 
respectively. Therefore, we have 
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E = E E E ( 3 - 27 ) 

PeS R+K {a 1 ,...,a R }c{l,2,...,R+K} P R es R p K es K 



Thus, we have the relation fl3.14|) , where a m 's correspond to j m 's 



Q.E.D. 
It is now clear that we obtain Theorem III.l from Proposition III. 2. 

IV. FORMAL BETHE STATE FROM THE ALGEBRAIC BETHE ANSATZ 

A. The algebraic Bethe ansatz for the XXX model under the P.B.C.s 

The formal Bethe state has been formulated in terms of the coordinate Bethe ansatz 
method in sec. 1. However, in this section, we show that it is also important in the context 
of the algebraic Bethe ansatz method. In fact, we show that the formal Bethe state \\M) 
with M rapidities corresponds to the state created by the B operators with the same set 
of rapidities. In the derivation, we use the method of the generalized two-site model first 
discussed by Izergin and Korepin |18|,|19 . 



Let us formulate some notation for the algebraic Bethe ansatz of the XXX model under 
the periodic boundary conditions. We define the L operator acting on the nth site of the 
one-dimensional lattice by 



L n (X) 



I XI n + Wn 2 V°n 



(4.1) 



y 2r/a+ XI n - T]a z n j 

Here we recall that the cr^'s are the Pauli matrices acting on the nth site. The monodromy 
matrix is defined by the product of L operators 

T(A) = L L {\ - gz)L x _i(A - q L ^) ■ ■ ■ L x {\ - q x ) (4.2) 

Here the free variables g^'s are called inhomogeneous parameters. Let us denote the operator- 
valued matrix elements of the monodromy matrix by 



T(A) 



1 ' A(X) B{\)^ 



V 



(4.3) 



/ 



C(X) D(X) 
The transfer matrix of the XXX model is given by the trace of the monodromy matrix 



r(X)=ti(T(X))=A{X)+D(X) 



(4.4) 



The Hamiltonian of the XXX model is derived from the logarithmic derivative of the ho- 
mogeneous transfer matrix where all the inhomogeneous parameters are given by zero. We 
note that when q^ = for k — 1, . . . , L, we call the transfer matrix homogeneous. Explicitly, 
we have 



tW~' A t(A ) 



1 



/ 1=1 



(4.5) 



\=r); qi=-=q L =0 ' 3- 

Let us consider the Yang-Baxter equations. With the R matrix 



/ 



R(X) 



1 
A 



A + 277 

2r] X 

A 2r/ 

A + 2r/ 

we can show the Yang-Baxter equations for the L operators 



V 



(4.6) 



R(X - n) (L n (X) ® L n {n)) = (LM ® L n {X)) R(X - n) 



(4.7) 



Here the symbol ® denotes the direct product of matrices. Applying eq. ( ^4.7| ) to each site, 
we can derive the Yang-Baxter equations for the monodromy matrix 



R{X - n) (T(A) ® T(fij) = (T(fj) ® T(A)) R(X - fi) 



(4.8) 



The Yang-Baxter equations for the monodromy matrices give the set of commutation rela- 
tions among the operators A(X), B(X), C(A), -D(A). For instance, we have 



A(X 1 )B(X 2 ) = f 12 B(X 2 )A{\ 1 ) - g 12 B(X 1 )A(X 2 ) 
D(X 1 )B(X 2 ) = f 2 iB(X 2 )C(X 1 ) - g 21 B(X 1 )D(X 2 ) 



(4.9) 
(4.10) 
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Here the symbols fu and gu denote the following 

Ai - A 2 - 2 V 2 V 

/l2 = r r , 912 = -T r (4.11) 

Al — A2 A\ — A2 

The operators A(X) and D(X) act on the vacuum |0) as 

A(A)|0) = a(A)|0), D(A)|0) = d(A)|0> (4.12) 

where a (A) and d(X) are given by 

a(X) = f[(X-q k + r ] ), d(X) = f[(X - q k - rf) (4.13) 

fc=i fe=i 

B. The Bethe ansatz wavefunction from the generalized two-site model 

Let us explicitly calculate the entries (or matrix elements) of the vector 
-B(Ai) • • -S(Am)|0). Let us consider two sub-lattices of the one-dimensional lattice with 
sites 1 to L: one consisting of sites 1 to n and the other of n + 1 to L. The monodromy 
matrix T(A) for the total lattice is given by the product of two monodromy matrices for the 
sub-lattices 

T(A)=T( 2 )(A)T«(A) (4.14) 

where T^(X) and T^(X) are given by 

T (1) (A) = L„(A-g n )L„_i(A - g n _i) • • -Li(A - qi) 

T^(A) = L L (X - q L )L L ^(X - q L ^) ■ ■ ■ L n+1 (X - q n+1 ) (4.15) 

Let us denote the elements of the monodromy matrices of the sub-chains by the following 



T (a) (A) 
Then, the product is given by 



( A (a \X) B^(X)^ 

yC^{X) D^(X) j 



for a = 1,2. (4.16) 
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A(X) B(X) 
C(X) D(X) 



\ ( 



AW(\) B® (A) 
C^(X) DW(\) 



\ 



( 



AW(A) fiW(A) 
CW(A) DW(A) 



\ 



(4.17) 



It gives four relations. Among them, we consider the following 



B(X) = A^(X)B^(X) + B^(X)D^(X) 



(4.18) 



By applying the formula ( f4.18| ) extensively, we can show the following [ ISpU 



M 



((2)) 



,(2) 



n € } ] (4.i9) 



n ^(A,)io) = e n b® n *s> io> n n /h n 

j=l 51,52 \keSi maS 2 J \keS 1 rneS 2 J \k&Si 

Here the symbol Y^Si L denotes the sum over all the combination of sets Si and 5*2 satisfying 
the condition: S1US2 = {1, . . . , M}. Here we have denoted B(Xk) by Bk, briefly. The symbol 
a k has been defined by the relation: A^ 2 \Xk)\0) — a k , and so on. We can prove the formula 
Q4.1SQ by induction on M. 

Let us consider the case where the total lattice of 1 to L is divided into r parts. Then, 
applying the formula ( 4.19Q , we can show the following |1~8|J20[ 



M ( (r')) / r 

UB(x,)\o)= e n n bw||o> 

j=l Si,...,S r \a=lm a £S ai 



x n 

l<a</3<r 



n n fkm 

\k&S a m£Sp 






,0») 



n ^ 



(4.20) 



Here the symbol Yl% ... s r denotes the sum over all the sets Si,...,S r which satisfy the 
condition: Si U 5*2 U ■ ■ ■ U S r = {1,2,..., M}. The formula ( |4.20| ) can be proven by induction 
on r. 

Let us consider the case of r = L where each of the sub-lattices is consisting of only one 
site. Then, we can show B k Bffl\0) = 0, for any k,m G S a . Thus, if any one of the sets 
Si, . . . , Sl contains more than one integer, then the contribution for the sets vanishes in eq. 
( |4.20| ). Therefore, we may consider only such partitions of {1, ... , M} into L sets: Si, . . . , Sl 
where each of S/s contains at most one integer. Let us take one such partition and assume 
that the symbols S X1 , . . . , S XM denote the non-empty sets with 1 < x\ < x 2 < • • • < Xm < L. 



21 



Then, there exists an element of the permutation group Sm such that S Xj = {Pj} for 
j — 1, . . . , M. Now we can derive the following expression from eq. fl4.20| ). 



M M 

U b(\ 3 )\o) = (zv) M E E n%lo) 

3=1 PeS M l<xi<-<x M <Lj=l 

M 

^ " ■ ; , 

,l<a</3<A/ / a=l \l<fc<j Q j a <k<L 



n /«ft» n n 4X n «s ( 4 - 21 ) 



Here we have used the following 

M M 

n4f|0) = (2r ? ) M n%|0) (4.22) 

3=1 3=1 

which is valid when each of the L sub-lattices consists of one site. Here we note that the 
product Ili<j<fc<M fpjPk is related to the amplitudes of the Bethe ansatz wavefunctions 

/ \ / \ _ \ j_ o \ H(P-' L 3-P- 1 k) 

n fpj*.= [ n fA n u-\-l ) (4 - 23) 

l<j<k<M \l<j<k<M j l<j<k<M \ A 3 Ak Z ''/ 

For the homogeneous case where qk = for k — 1, . . . , L, we have 

n ( n 4i n .ffi) = n {^^ (r^fl (424) 

Putting 1] = —i, we have 

M A'/ / \ / M 

ni?(A J )io)=F 1 ({A J }) e e rKi°> n w- e XP u^^ 

j=l Pe5 M l<xi<-<a:M<ii=l \l<a</3<M / \ j=l 

(4.25) 



where the factor Fi({Xj}) has been given by 



^({A,-}) = (2rj) M g %-^T (4-26) 

7=1 A i "T 2 



Noting the relation (|4.23| ) we obtain 



M 

Y[B(\ J )\0) = F 1 ({\ J })F 2 ({\ J }) 
j=i 

(M \ / Af 

I] er" |0) A M (P) [Ai, • ■ ■ , Am] exp i E k Pj Xj ] (4.27) 
3=1 ) \ 3=1 
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where the factor F 2 has been given by F 2 ({Xj}) = Tii<j<k<M fjk- hi terms of the notation 
of the formal Bethe state \\M), we have 

B(X l )...B(X M )\0) = \\M) x F 1 ({X j })F 2 ({X j }) (4.28) 

Thus, we have shown that the formal Bethe state \\M) is nothing but the vector 
-B(Ai) • • ■ B(Xm)\0) where all the rapidities are given by free parameters. Here we remark 
that the derivation of ( |4.21| ) for the cases of M = 1 and M = 2 has already been given 



explicitly such as in Refs. [1,10] ■ 

We give some remarks. The expression ( |4.21|) also holds for the XXZ model, where fjk, 
a(A) and d(X) are replaced by 

sinh(Aj — A fc — 2rf) sinh(2r]) 



sinh(Aj — Afc) sinh(Aj — A&) 

L L 

o(A) = J] sinh(A -q k + rj), d{X) = H sinh(A -q k -rj). (4.29) 

fc=i fc=i 

Corresponding to eq. ( |4.23| ) , for the XXZ model we have 

n fpjpk=[ n f») n (¥) ( 4 - 3 °) 

l<j<k<M \l<j<k<M I l<j<k<AI \JJ k / 



A similar relation with (|4.21| ) has also been derived for the algebraic Bethe ansatz of the 



elliptic quantum group [21], where fjk are expressed in terms of the elliptic theta functions. 



V. THE SPECTRAL FLOW UNDER THE TWISTED B.C.S 

A. Formal Bethe states under P.B.C.s 

Let us discuss the formal Bethe state from the viewpoint of the ground-state solution 
of the XXZ model given by Yang- Yang M . We consider the case of the anti-ferromagnetic 
Heisenberg model. In this section we assume that L is even. We note that the relation 
between momentum and rapidity is slightly different from that of seel and sec. 2 due to 
the gauge transformation p. 
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The Hamiltonian of the XXZ model is given by 

« = " J E {SfSf +1 + Sj[SJl +1 + A5?5? +1 ) . (5.1) 

where A is called the anisotropy parameter. Let us discuss the anti-ferromagnetic XXX 
model. Hereafter we assume A = — 1. Under the periodic boundary condition: S^+i = Si, 
we have the Bethe ansatz equations 

exp(^) = (-l)^ n eX P^^ + ^] + l-2(-l)exp(,A;,) 



\ m exp[i(kj + h)] + 1 - 2(-l) exp(ifc^) ' 
for j = 1, . . . , M. (5.2) 

We consider the case where all the momenta k/s are real. Taking the logarithm of the Bethe 
ansatz eqs. ( |5.2|) , we have 

M 

Lp J =2nI J - Yl ®(Pi,Pe), for j = l,...,M. (5.3) 

Here Ij = (M — l)/2 (mod 1). The function 0(p, q) has been given by || 

Q(v a) ~ 2 tan- 1 ( (-l)sin((p-g)/2) \ 

0(p,g)-2tan { cos{{p + q)/2) _ { _ l} cos{(p _ q)/2) ) (5-4) 

In Ref. ||, some important analytic properties of the function Q(p,q) have been given. In 
particular, we notice the following: When — ir < p < it and — 7r < q < ir, we have 

lim Q(tt — e, q) = — tt, lim Of— 7r + e, o) = 7r, (5.5) 

e(p,q) = -e(q,p). (5.6) 

Here we note that e should be positive in the limits. In Ref. ||, the function is mainly 
discussed within the range (— n, n) for p and q. 

Let us first show that the limit: p — > —it under p > —it corresponds to the limit of infinite 
rapidity: A — > oo discussed in the previous sections. Let us introduce the transformation 
between momentum k and rapidity v for the anti- ferromagnetic case (A = —1) 

exp(*fc) = (-1) (^) . (5.7) 
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Here the extra (—1) factor corresponds to the gauge transformation, which is consistent 
with the periodic boundary condition if L is even. By taking the logarithm of the relation 
Q5.7| ) and choosing the branch of the logarithmic function such that is consistent with the 
function 6(p, q), we have 

k = -2 tan" 1 *;. (5.8) 

It is clear that momentum k is in the range (— tt, tt) if v is finite. If we take the limit v — > oo, 
then the momentum k approaches to — tt from the above (k > —tt). Thus, we can extend 
the range of real momenta as follows 

- tt < k < tt. (5.9) 

Let us call momentum kj regular if it satisfies the condition: — tt < kj < tt. We now 
introduce a symbol 7\ for (L — M — l)/2. From the viewpoint of the string-hypothesis, it is 
shown that the integer (or half-odd integer) Ij for a regular real momentum kj satisfies the 
condition: \Ij\ < T x [0,|||13],[22 . 



We now consider such a solution to the Bethe ansatz equations that has the momentum 
k Q = —tt. Suppose that the set of R regular momenta k\ , . . . , kpt and one non-regular 
momentum k = —tt gives a solution to the Bethe ansatz equations ( |5.3|) with M = R + 1 
down-spins. We have the following equations. 

R 

Lkj = 2TTlj- Y, ®(kj,k e ), for j = 1,...,R, (5.10) 

R 

Lfc = 27r/o-£e(fc ,£;,). ( 5 - n ) 

Here we have defined the equation for the momentum ko = —tt by taking the limit of k to 
— tt under the condition k > —tt. We also assume that Ij = R/2 (mod 1) and \Ij\ < T\ — 
(L - R)/2 - 1 for j = 1, . . . , R. (M - 1 = R.) 

Let us construct solutions to eqs. ( |5.10| ) and ( |5.11| ). Making use of the limit ( p.5| ), we 



can easily show that ko = —tt gives a solution to the equation ( J5.ll ) with I = — T\ — 1, 



where T\ = (L — R)/2 — 1. Recall that for ko = —tt, we assume the limit: ko — > —tt under the 
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condition: k$ > —it. From the property (|5.5|) , we can show that the eq. ( |5.10|) is equivalent 

to 

1 R 

Lk j = 2n{I j + -)- Y. ®(kj,k e ), for j = 1, . . . ,R. (5.12) 

1 e=i,i& 



We note that the set of equations ( |5.12j ) is equivalent to the standard Bethe ansatz equations 



Q5.3D for R down-spins, where their numbers (/.,■ + 1/2) 's satisfy the parity with R down-spins: 
17 + 1/2= {R-l)/2 (mod 1). 

Let us summarize the construction of solutions to (|5.10|) and ( |5.11|) . We first consider the 
standard Bethe ansatz equations ( |5.3|) (or (|5.12|) ) with R down-spins where their quantum 
numbers J/s are given by Ij = Ij + 1/2 for j = 1,...,R. We denote the solutions by 
fci, . . . , kn for Ii, . . . , I R , respectively. Then, the solutions to eqs. ( |5.10| ) and (|5.11| ) are given 
by 

kj = kj, for j = 1, . . . , R, 

k = —ii, (5.13) 

where J = -{L-R)/2. 

It is interesting to note that the set of momenta ko, . . . , k R , constructed in the above 
corresponds to such a solution of the Bethe ansatz equations that gives the non-regular 
eigenstate \R, 1). In fact, we obtain the non-regular eigenstate \R, 1), by taking the limit of 
the formal Bethe state | \R, 1; A) with the momenta ko, . . . , kji, where k corresponds to the 
infinite rapidity. 

B. Bethe states under the twisted B.C.s with a small twist 

Let us now consider the twisted boundary conditions 

^ +1 = 5 1 ± exp(±z$), S z L+l = S{. (5.14) 

Here we call the variable $ the twisting parameter. The Bethe ansatz equations under the 
twisted boundary conditions are given by p3| - 



26 



exp(zL^) - ( 1) «p(,*) J^ ^.^ + ^ + l _ ^ ^.^ , 

for j = l,...,M. (5.15) 



Let us assume that all the momenta pj's are real. Taking the logarithms of eqs. (|5.15|) , we 
have 

M 

Lp J = 2TrI j + $- Yl ®(PhVl), for j = l,...,M. (5.16) 

where I j = (M - l)/2 (mod 1). 

Let us discuss the adiabatic behaviors of solutions of the Bethe ansatz equations ( |5.16| ) 



under the twisted boundary conditions, where $ is a very small positive number, and is 
increased adiabatically. Here we assume that R + 1 momenta: po,Pi, ■ ■ ■ ,Pr are solutions 
of the Bethe ansatz equations where one momentum p is close to —n (p ~ — 7r). We also 
assume that R < L/2. We consider the following equations. 

R 

Lp j = 2ixl j + ®- Y, e (P^), forj = l,...,i?, (5.17) 

R 

Lp = 2nl + $ - Y,Q(po,p t ). (5.18) 

e=i 

Here, Ij = R/2 (mod 1) and \Ij\ <T X = (L- R)/2 - 1 for j = 1, . . . , R. (Note M - 1 = R.) 
Now, let us introduce a small positive number e and express po as po = — vr + e. Then, we 
can expand the function Q(p, q) in the following 

0(-7T + e, q) = n - 2e + e 2 tan - + 0(e 3 ), 
0(7r - e , q) = -7T + 2e + e 2 tan | + 0(e 3 ). (5.19) 



Substituting the expansion into eq. ( |5.18| ) we have 



27t/ + $ = -tt(L -R) + e(L - 2R) + 0(e 2 ). (5.20) 

Let us consider the adiabatic change of the twisting parameter $. When we change the 
parameter $ infinitesimally, then the quantum number Iq does not change at all. Thus, we 
obtain the following solutions 
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I = -L- R/2 = -Ti - 1, 
e = $/(L-2 J R) + 0(e 2 ). (5.21) 

The solution: p = —n + $>/(L — 2R) + 0(e 2 ) can be considered as a regular solution, since 
it satisfies the condition —n < po < —tt. Furthermore, it is clear that po approaches to 
fco under the limit $ — > with $ > 0. Therefore, we conclude that the formal solution 
ko = —ix under the periodic boundary conditions corresponds to the regular solution po = 
— ti + $/(L — 2R) + 0(e 2 ) under the twisted boundary conditions. 

We make a remark on the number of solutions to eq. ( |5.16 ). Under the periodic boundary 



conditions, there are 27\ + 1 = L — M regular solutions. Under the twisted boundary 
conditions, however, we have one more regular solution and 2Ti + 2 = L — M+l regular 
solutions in total. 

C. Analytic continuation of &(p,q) and 27r-shift of momentum 

In the previous subsection we have considered the case where the twisting parameter $ 
is very small. Hereafter, we discuss how the ground-state solution changes with respect to 
the parameter $, and then we shall show that the spectral flow should have the period of 
47r, at the end of sec. 5. 

Let us consider an extension of the function Q{p, q) of the variable p defined on the range 
(—it, it) into a continuous function defined over (— oo, oo). Taking into account the analytic 
property (|5.5| ), we make the extension with respect to p as follows 



E(p, q) 



-(2n+l)7r, if p = (2n + 1)tt, foraninteger n, 



e(p - 2tt 



P + TC 



2vr 



,g)-27r 



P + TT 



2tx 



otherwise. 



(5.22) 



Here the symbol [x] denotes the Gauss' symbol. We recall that p is an arbitrary real number 
in eq. (|5.22|) . We also extend the function 0(p, q) with respect to q by assuming the relation 

o(p,q) = -o(q,p)- 

In terms of the extended function, the Bethe ansatz equations are given by 
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M 

Lp j = 2nI j + $- ]T S 0w), f o r J = l,-.-,M, (5.23) 

where Ij = (M — l)/2 (mod 1). We recall that the range of pj's are given by the extended 
zone (—00, 00) in eqs. ( |5.23| ). We note that for regular solutions, eqs. ( p.23| ) are equivalent 



to the standard Bethe ansatz equations ( |5.16| ). We recall that when pj's are regular, then 
they satisfy — it < pj < 71. 

Let us discuss the number of all the possible solutions to eqs. ( |5.23|) . We assume that 
Pi, . . . ,Pm are solutions to ( |5.23|) with the quantum numbers Ji, . . . , Im, respectively. Let 



us take a suffix j\ from 1, . . . , M. We consider momenta pj given by 

p n = p h + 2vr and p 3 = Pj if j ^ 3l (1 < j < M). (5.24) 

Putting pj 's into the eq. ( |5.23|) , we see that they give solutions to eq. ( |5.23| ) with the 



quantum numbers Ij's where they are given by 

I jl = Ij.+L-M + l, and Ij = Ij + 1 for j =£ ji{l < j < M). (5.25) 

Thus, we may regard the number L—M+l as the period of the quantum numbers Ij's. There- 
fore, we may consider only L — M + l different values for J^-'s such as —T\ — 1, — 1\, • • • , T\. 
This choice is consistent with the range of momentum — tc < k < ir. (Iq — —T\ — 1 corre- 
sponds to k = —it or v = 00 when $ = 0.) Thus, if there is a one-to-one correspondence 
between kj's and Ij's, then the number of solutions of the Bethe ansatz equations is given 
by L — M + 1. We note L — M + 1 = (2Ti + 1) + 1, where 27\ + 1 is the number of regular 
solutions given by the string hypothesis. 

Let us make a remark on the change in the total momentum P to t induced by the shift: 
Pj 1 — > pj 1 + 2tt. In fact, it is consistent with the shift of the quantum numbers Ij's. We 
denote by Ptot the total momentum for the momenta pj's and Ij's. Then, we have 

Ptot = Y,Pj = 27r + &j; 



3 3 



2tt ^ f 2tt _ I-+L-M + 1 n 2vr _ 

X^ J i = T ^ (/j " 1} + Z = 2vr + T^ Jj - 



29 



(5.26) 



D. Period Air of spectral flow under the twisted B.C.s 

Let us discuss how the ground-state solution changes when we increase the twisting 
parameter $ adiabatically from to Atc. In fact, it is known that the spectral flow of the 
ground state energy has the period of 4ii with respect to the twisting parameter [p5|-|29]. In 



the following paragraphs, we shall show that the shift ( |5.25 ) on the I/s is consistent with 



the 47r period of the spectral flow of the ground-state energy. 

We consider the half-filling case where M = L/2. Here, we have 7\ = (L/2 — l)/2. Let 
symbols Pj{<&) for j = 1, . . . , M denote the momenta satisfying the Bethe ansatz equations 
( f5.23| ) under the twisted boundary condition. We assume that when $ = 0, the momenta 
are given by those of the ground state, where —n < Pj{0) < tt and Ij = —T\ + (j — 1) for 
j = 1, . . . , L/2. We may assume that pi(0) < ■ • • < Pl/2(0). 

Let us recall the adiabatic hypothesis that when we change the parameter $ infinitesi- 
mally, then the quantum number Ij for momentum Pj($) does not change at all, while the 
momentum Pj(&) changes proportionally to the infinitesimal change of $. We note that 
under the adiabatic hypothesis we have fci($) < • ■ • < ki/2^)- Now we consider the case 
when $ is very close to 2ir. We set $ = 2n — 5, where 5 is a small positive number. Making 
use the expansion ( |5.19|) , we can solve the Bethe ansatz equations (|5.23| ), and we obtain the 
following 

Pl/2 (2ti-5) =ti- 5/2 + (5 2 ). (5.27) 

Thus, when $ = 2tt, we have momenta pi(2ir), . . . ,pi/2(2vr), where their quantum numbers 
are increased by 1 and they are given by 

J i = -Ti + 0'-l) + l for j = I,..., L/2. (5.28) 

Here we note that pl/ 2 (2ii) = -a and Il/2 = Ti + 1. Let us now apply the shift ( |5.25| ) to the 
system of the solutions. We note that pl/2{27t) — 2tt = — it and Ti+1 — (L — M+l) = — Ti — 1. 
Thus, they are equivalent to the following set of solutions to (p.23[) 



Pl/2(2tt) - 2tt = -tt, Pi(2tt), . . . ,p (L _ 2 )/2(27r) 
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with their quantum numbers given by 

h/2 = —T\ — 1, h = — Ti, . . . , Ti — 1, 

respectively. We further increase the twist parameter $ until it becomes Air. Then, we see 
that the set of momenta 

Pl/2(4tt) - 2vr, pi(47r), . • • ,p L /2-i(4tt) 

satisfy the Bethe ansatz equations ( |5.23|) with 



h/2 = —Ti, h = —T\ + 1, • • ■ , I(l-2)/2 = T\. 
Thus, we obtain 

PL/2(4vr) - 2VT = Pi(0),Pi(47r) = p 2 (0), • • • , P(L-2)/2 (47r) = P(L~2)/2- 

The set of momenta for $ = 4ir is equivalent to that of $ = 0. Therefore we conclude that 
the solutions to the Bethe ansatz equations which corresponds to the ground-state at $ = 
has the period of 4n. 

VI. DISCUSSIONS 

In this paper, we have explicitly shown that any non-regular eigenvector is derived from 
the Bethe ansatz wavefunction with infinite rapidities, for the one-dimensional XXX model 
under the periodic boundary conditions. The formula (|2.10| ) for the amplitudes of the Bethe 
ansatz wavefunction has played a central role in the proof. 

Let us discuss the string hypothesis, explicitly. It is based on the assumption that the 
Bethe ansatz equations (|2.8|) have complex solutions given in the following 

v n J = v n a + i(n + l- 2]) + e n J for j = 1, . . . , n. (6.1) 

Here, it is also assumed that the absolute values of the correction terms |e" J | should be 
very small. The set of complex rapidities v^ for j = 1, . . . , n is called an n-string solution 
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PI JTSlJT^I . The value v™ is called the center of the string solution. The number n is called 



the length of the string solution. 

Let us discuss the formula ( |2.10| ) of the amplitudes from the viewpoint of the string 



hypothesis. For any given n-string solution, we set the n rapidities in the string in such 
an order that v £ — v* « 2i(k — j) for any j < k with 1 < j, k < n. Then, the value of 
the amplitude Am{P) given by eq. Q2.10D becomes stabilized and well-defined, since we can 
avoid the appearance of any very small factor of 0(e) in the denominator of eq. ( p,10|) . 



Let us discuss the number of solutions to the Bethe ansatz equations for the strings of 
length n []I|,|l3f . Let us define number T n by 

T n = l[N-l-jr i t nm M m ) (6.2) 

Z \ m=l I 

Here M m denotes the number of string solutions of length m and t nm is given by 

t nm = 2min(n, m) - 5 nm . 

Under the periodic boundary conditions ($ = 0), it is discussed that the number of string 
solutions of length n is given by 2T n + 1 |T3||[j| . We can show that if there are 2T n + 2 
different string solutions of length n, then all the solutions to the Bethe ansatz equations 
correspond to a complete set of the eigenvectors of the XXX Hamiltonian under the twisted 
boundary conditions. The result of the present paper suggests that the K infinite rapidities 
of a non-regular eigenvector \R, K) might correspond to a i^-string solution under the twisted 
boundary conditions. Thus, we have a conjecture that any non-regular eigenvector under 
the P.B.C.s of the form \R, K) should correspond to a regular eigenvector with a iT-string 
solution under the twisted B.C.s. It seems that the conjecture should be consistent with 
the result of Ref. [|I7]]. However, a detailed numerical research on A'-strings with large K's 
should be performed such as studied in Ref. |30| . 



Finally, we give a remark on a possible application of the result of the present paper to 
the XXZ and XYZ models. Recently, it has been shown that under the periodic boundary 
conditions, the one-dimensional XXZ Hamiltonian at the q root of unity conditions has the 
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sli loop algebra symmetry [pl|-|33||. In fact, we can discuss the spectral degeneracy of the 



XXZ model at the root of unity conditions in terms of the algebraic Bethe ansatz method 



by applying some of the techniques developed in the paper |34| : combining the expression 
( |4.21| ) with the formula ( |4.30| ) of the amplitudes Am(P)'s, we can construct singular solutions 
related to the s/2 loop algebra. Thus, we can show the validity of the construction of the 



complete A-string solutions discussed in Ref. |33| in the level of eigenvectors. We can also 
prove it by showing that the limits of the Bethe ansatz wavefunctions satisfy the sufficient 
conditions for the eigenvectors of the XXZ model, which are summarized in Appendix E. 
Surprisingly, a similar method can also be applied to the analysis of the spectral degeneracy 



of the XYZ model addressed in Ref. |31|. The details will be discussed in subsequent papers. 
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A. APPENDIX: FORMULA FOR THE ACTION OF SPIN-LOWERING 

OPERATOR 

Let us introduce some symbols. First, we shall abbreviate the symbol J2i<xi<-<x m <l 
by J2xi<-<x M i * n sn °rt- Second, for a non-negative integer K, we denote by the symbol 
^<{ji,j2,-,jm}c{i,...,m+k} the summation over all the subsets {ji, j 2 , ■ ■ ■ ,Jm) of {1, 2, . . . , M + 
K}, where jVs are set in increasing order ji < ■ ■ ■ < Jm- Thus, the two symbols in the 

33 



following express the same sum. 



E = E (A-i) 

{ji,...,j M }c{l,...,M+K} l<ji<-<j M <M+K 



Proposition A.l Recall that \M) denotes an arbitrary vector with M down-spins defined 
by eq. t \2.3j ). We denote by \M,K) the vector obtained from \M) multiplied by the power of 
the spin-lowering operator 

\M,K) = ^(S tot ) K \M) (A.2) 

Then, we can show the following formula 

\M,K)= J2 ( E 9(x jl ,...,x JM ))cr- 1 a^---cr- M+K \0) (A.3) 

xi<-<x M+K \{ji,...J M }c{l,...,M+K} J 

(Proof) We prove the formula ( |A.3j ) by induction on K. 

(i) We show ( |A.3| ) for the case of K — 1. Applying S^, t to |M), we have 



L 

S tot\ M ) = E a <7 E 9(xi, ■■■, x M )a^ ■ ■ ■ cr-JO) 

y=l Xi<--<Xm 

~ (y<x\ y<x 2 L \ 

= E E + E+---+E ) 9{xi,...,x M )or;o:- i ---a~ M \0). 

xi_<—<xm \j/=i y> x i v>xm) 

(A.4) 

We note the following calculation. 

~ y<xj+i 
E E 9(xu ■■■, x M )(y y o- x ■ ■ ■ a~ M 

xi<---<x M y>Xj 

E 9(%i, • • • , %)^ ■ ■ ■ 0-<?-<T- j+1 

x 1 <--<Xj<y<Xj + i<--x M 



j"y " x j+l X 'M 



E 9{xi, ■■■, x j7 x j+2 , . . . , x M +i)o- xl ■ ■ ■ a (A.5) 



Xl< — <XM+1 



(j+l)th,...,Mth 

In the last line, we have replaced the symbols y, Xj+i, . . . , and xm by Xj+i, £j+2, • • • , and 
Xm+i, respectively. Substituting (|A.5| ) into ( |A.4| ), we have 
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s tot\ M ) = E {g{x 2 ,...,XM+i)+g{xi,X3,...,x M +i)^ 

xi<---<a; A / + i 

+ • • • + g(xx,x 2 , ..., x M ))o-~ • • • (t~ m+1 |0) 
= E E 9(x h , ..., x jM ) U" • • • <r~ f+i |0). 

x 1 <-<x M +i \{iivjAf}C{l,2,...,M+l} / 

(A.6) 



Thus, we have the expression (|A.3| ) for the case of K = 1. 

(ii) Let us assume the expression ( |A.3j ) for the case of K. Then, we show the case of K + 1 

in the following. 

S M \M,K) 

= J2 a y[ E E 9(xh ,■■■, Xju)°Zx ■ ■ ■ a M+K |0) 

y=l \xi<-<x M+K {j lt ...,j M }c{l,...,M+K} ) 

~ (y<x\ y<X2 L 

= E E + E+---+ E 

xi <—<x m+k \y=l y>x\ y>x M +Kj 

x E 9(x jl ,...,x JM )a y a- i ---a- M+K \0). (A.7) 

{Ji-Jm}c{1,...,M+K} 

By a similar method for the case (i), we can show the following 

~ y<x l+1 

E E E 9{xn,--.,x JM ) (j-a~ x ■ ■ ■ a~ M 

X!<-<x M+K y>x t {j lt ...J M }c{l,...,M+K} 



E E 9 {x h ,..., x JM ) a xi ■ ■ ■ o Xi a y a X(+i 

x 1 <-<x l <y<x l+1 <--<x M +K {j 1 ,...,j M }c{l,...,M+K} 



°x M 



E E 9{x h ,...,x JM )a Xi ---a XM+i (A.8) 

x\<-<x M+K +i {ji,---,j M }c{i,...,e,e+2,...,M+K+i} 

In the last line, we have replaced the symbol y by Xt+i, and £+1, . . . , M by £ + 2, . . . , M + l, 

respectively. Substituting (|A.8| ) into ( |A.7| ), we have 



E f E + E +•••+ E 

xi<-<x M +K+l \{j 1 ,-j M }c{2,3,...,M+K+l} {ii,--j M }c{l,3,...,M+K+l} {j 1 ,-,j M }c{l,2,...,M+K}; 



xg(x h , . . . ,x JM )a xl - ■ -a XM+K+i \0) 



(K + l) E E 9(x jlJ ... 7 x JM ))a- l ---a- M+i \0) (A.9) 

xi<-<x M +K+l \{j 1 ,...,j M }c{l,2,...,M+K+l} 
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In the derivation of the last line, we note that after selecting M integers ji,J2, ■ ■ ■ ,Jm from 
the set {1, 2, . . . , M + K + 1}, there are (K + 1) ways for choosing one more element from 
the remaining K + 1 integers. Thus, we have the factor (K + 1). 

Q.E.D. 

B. APPENDIX: FORMAL BETHE STATE WITH THREE INFINITE 

RAPIDITIES 

Let us discuss the infinite limit of the formal Bethe state \\R, K; A) of the case R = 
and K = 3. Here, Vi, v 2 and v 3 are additional rapidities given by v± — A + Si, v 2 = A + S 2 
and t>3 = A + ^3. Let us denote Si — S 2 and S2 — S3 by A12 and A23, respectively. After taking 
the limit of sending A to infinity, we have 

^123(00) = 1, ^132(00) = — -, A 2 i 3 (00) 



423l(°o) 

^312(oo) 



A 23 + 2i ' -'- " ' A 12 + 2i' 
A 12 -2i\ /A 12 + A 23 -2^ 



Ai2 + 2i/ VAi2 + A 23 + 2i 
A 12 + A 23 -2i\ /A 23 -2i 



Ai2 + A 23 + 2i7 VA 23 + 2i/ ' 



a ( \ ^ A 12 - 2i \ f A u + A 23 - 2i \ ^ A 23 - 2i \ 

A32l(oo) = taaJ U2+A23+2J toaJ • (R1) 

The limit of the Bethe ansatz wavefunction for the formal Bethe state 1 10, 3; A) is given by 

lim f$(xi,x 2 ,x 3 ;ki(A),k2(A),k 3 (A))=C 3 (B.2) 

where the constant C 3 is given by 

C 3 = ]T A 3 (P)[Si,S 2 ,S 3 ] 
PGS3 
= x A 23 - 2i A12 - % fA u - 2i\ ( A12 + A 23 - % 



+ 



A 23 + 2i A12 + 2z VA 12 + 2i) V A 12 + A 23 + % 
A12 + A 23 - 2i\ /A 23 - 2i\ /A12 - 2i\ /A12 + A 23 - 2i\ /A 23 - 2i 



A12 + A 23 + 2i/ V A 23 + 2i/ V A12 + 2i) V A12 + A 23 + 2i/ VA 23 + 2% 



(B.3) 



Thus, we have 
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A lim 1 10, 3; A) = C 3 £ a~ t a~ 2 cr" |0) 

l<Xl<X2<X3<L 



Cs^S^flO) 



3! 
C 3 |0,3) (B.4) 



C. APPENDIX: SOME USEFUL PROPERTIES OF THE SYMMETRIC GROUP 



We introduce some notation of the symmetric group [BH] . Let M be a positive integer. 



We consider the permutation group Sm of integers 1,2,..., M. Take an element P of Sm- 
We denote the action of P on j by Pj for j = 1, . . . , M. Let us introduce the symbol 
(ziz 2 • • • v) of the cyclic permutation where ij is sent to z,- + i for j = 1, ... ,r — 1, and i r is 
sent to i\. It is known |3(|] that any permutation P can be decomposed into a product of 
disjoint cycles such as follows 

P = (*ii 2 ■ ■ • v)(jij 2 • ■ ■ js) • ■ ■ (^i4 • • -4)- (C.i) 

Here, any two of the cycles share no letter (or integer) in common. The factorization (|C . 1|) 



is unique except for order of the factors [36 



For a given permutation P with a factorization of disjoint cycles such as eq. ( CJ.1 ), we 



denote by N(P) the sum (r — 1) + (s — 1) + ■ ■ • + (u — 1). Then, we can show that the parity 
of the permutation P is equal to that of N(P). Hereafter, we shall write by the symbol 
a = b (mod 2) that integers a and b have the same parity. We first recall that the cycle 
[i\i2 ■ ■ -ir) can be written as the product of r — 1 transposition such as 

(hi 2 ■ --ir) = (iii r ){hir-i) ■ ■ ■ (hk)- 
Thus, the parity of the cycle is given by that of r — 1. Let us denote by the symbol e(P) 



the sign of permutation P. Then, we have [pq] 



; (p) = (_i)(r-i)+(s-i)+-+0-i) _ f_i\N(P) (Q_2) 
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Let us introduce ordered pairs of integers. We take two different integers j and k, and 
consider an ordered pair < j,k >. We distinguish < j,k > from < k, j >. Let us consider 
the action of a permutation on ordered pairs. We take a permutation P of Sm and two 
integers j and k satisfying 1 < j < k < M. We denote by < Pj, Pk > the action of P on 
the pair < j,k >. If Pj > Pk, we call the pair < j,k > is transposed by P. 

Let the symbol T(P) denote the number of all such pairs < j,k > that are transposed 
by P among all the ordered pairs < j,k > with the condition 1 < j < k < M. Then, we 
can show the following. 

Lemma C.l The parity of an element P of Sm is equivalent to that of the number T(P): 

N(P)=T(P) (mod 2). (C.3) 

(Proof) We now prove the lemma based on induction on M of Sm- It is easy to see that when 
M = 2 the statement is true. Let us now assume that eq. ( |C.3| ) holds for all permutations 
P of Sr if R < M. Let us take an element P of Sm- Then, we may assume that the 
permutation P has a factorization of disjoint cycles such as shown in eq. ( |C.1| ). Suppose 
that P has the same factorization with eq. ( |C.1| ). We take a cycle (i\i2 • • • i r ), which is one 
of the disjoint cycles, and we denote by B the set {ii,i 2 , . . . ,i r }- We also denote by Sm the 
set of M integers: T, M = {1, 2, ... , M}. We now consider the subset A of the set S A/ that 
is complementary to the set B: A = Hm — B. We define permutation Pa by 

PA=(jlJ2---js)---(£lt2---Q- (C.4) 

Note that Pa is a permutation of A and it does not change any letter in B: PaIj = ij for 
j = 1, . . . ,r. Thus, we see that T(P) and T(Pa) + (r — 1) have the same parity. Here, we 
note that T((i\i2 ■ ■ -i r )) = r — 1, i.e., r — 1 pairs of the elements in B are transposed by 
(ii«2 • • 'ir), an d therefore by P. On the other hand, since Pa is a permutation of A, it is 
equivalent to an element of Sm-t- From the induction hypothesis, we have that N(Pa) and 
T(Pa) have the same parity. Thus, we have 
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T(P) = T{P A ) + (r - 1) (mod2) 
= N(P A ) + (r - 1) (mod2) 

= N(P) 

Therefore, T(P) and N(P) have the same parity. 

Q.E.D. 
We now have the following . 
Proposition C.l Let P be an element of Sm- Then, we have the following identity. 

e (P)= H (-i^p-^-p- 1 *) (C.5) 

l<j<k<M 

(Proof) Let us note the following 

T(P)= J2 H{p- l j-P- l k) (C.6) 

l<j<k<M 



Then, we can show eq. ( |(J.5| ) from the previous lemma and eq. ( C.2 ). 



D. APPENDIX: PROOF OF THE "PAULI PRINCIPLE" 

We give a simple proof for the "Pauli principle" of the Bethe ansatz that when there 

are two rapidities of the same value, then the Bethe ansatz wavefunction of the XXX model 

vanishes. We note that it is also proven by the algebraic Bethe ansatz method in Ref. 

2Cfl . However, the proof in this appendix is much more elementary; it is only based on the 



expression ( |2.10| ) of the amplitudes A M (P)'s. In this appendix, we assume that rapidities 



vi . . . , vm are free parameters. 

Let us take a pair of integers a and b such that 1 < a < b < M. Then, we show that the 
Bethe ansatz wavefunction f M with the amplitudes defined by eq. ( |1.3| ) (equivalently by 
eq. ( |2.1(J| )) vanishes if k a = k b (i.e., v a = Vb). Let the symbol (ab) denote the permutation 



between a and b. Then, we have 
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M 

fw (xi, ■ ■ ■ , x M ] kx,..., k M ) = Yl a m(P) exp(i ^ k Pj Xj) 

PeS M 3=1 

= 9 Z A M (P) exp li f] kpjXj + - Z A((ab)P) ex P I » J2 h {oh)P)j Xj (D.l) 

z Pe5 M V J =1 / Pg5 m V i=i / 

Here we have replaced P by (ab)P in the second term. Considering the cases when j = P -1 a 

and j = P~ l b, we can show 

M M 

Y kpjXj = k a x P ~i a + k b x P ~i b + Y kpjXj (D.2) 

3=1 j=X;j=£P- 1 a,P- 1 b 

M M 

/ , k(( a f))p)jXj = k( a b) a Xp-i a + k( ab ) b Xp-i b + 2_^ k( ab )pjXj 

3=1 j=X;j^P- 1 a,P- 1 b 

M 

= k b x P -i a + k a x P -i b + Y k Pi x J ( D - 3 ) 

j = X;j^P- 1 a,P- 1 b 

When k a = k b = k, we have 

fif\x 1 ,...,x M -k 1 ,...,k M ) = l Y (A M (P)+A M ({ab)P))) 

1 PeS M 

I M \ 

x exp ik(xp-i a + xp-i b ) + i J^ kpjXj (D.4) 

\ j=X;j^P- 1 a,P-n J 

We now show that Am(P) + Am ((aft) P) = for any P G Sm- Here we introduce the 
following symbols 

e(j, k) = V i- Vk ~ 2 * H{j, k- P) = HiP-'j - p-'k) (D.5) 



Then, the amplitude Am(P) given by eq. Q2.10D is expressed as 



A M (P)= II e(j,k) H W>V (D.6) 

!<3<k<M 

Let us consider the six cases for the integers j and k in the above product: j = a and k = b; 
j < a and k = a; j < b and k = b where j ^ a; j = b and k > b; j = a and k > a where 
k 7^ b; j 7^ a and k ^b. We have the following 

A M {P) = e(a, bf^ J] e(j, a) H ^ g e(j, &)*&V> 
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M M 

x J] e(b,k) H ^ k ^ n e(a,k) H ^ k ^ ]J e(j,k) H ^ k ^ 

k=b+l k=a+l;k^b l<j<k<M;j,k=£a,b 

e(a,a) H ^ b ^l[e(j,a) H ^ a ^ +H ^ a ^ [J e(j,a) H ^ 

j=l j=a+l 



M 6-1 

x A e(a,j) P(6j;P)+P(aj;P) II e(a,j) H(aj;P) II <h k) H ^ k ^ 

j=b+l j=a+l l<j<k<M;j,k=£a,b 

= (_1)H(°.M>) J] e(a,j)^ aj;P) - P(j ' 6;P) x JJ e(j, fc)^.M>) 

j=a+l l<j<k<M;j,k^a,b 

a-1 M 

x J] eO» ff0>;P)+ ^' 6;P) II e(o,j) P(6j;P)+P(a - ?;P) (D.7) 

j=l j=6+l 

Here we have used the relations e(j,a) = e(j,b), e(a,j) = l/e(j,a), e(a,b) = e(a,a) = — 1, 
and so on. In a similar way, we have 

A M ((ab)P) = (-l) H ( b <« p ) H e(a,j) P(6j;P) - P(j ' a;P) x [] e(j, A;)^ fc;P ) 

j=a+l l<j<k<M;j,kj^a,b 

a-1 M 

xne(j,< (wP)+W) II e ( a>J -)«(W^)+H(aj;P) (D . 8) 

3=1 3=b+l 

Noting the relation: H(j, k; P) — 1/2 = —(H(k,j; P) — 1/2), we can show 

H(a, j- P) - H(j, b; P) = H(p- l a - P^j) - H{p- l 3 - p- l b) 

= -H(-P- l a + P- 1 ]) + H(-P- l j + P- X b) 

= -H(j,a;P) + H(b,j;P). (D.9) 

Thus, we have 

A M (P) + A M ((ab)P) = ((-l) H ^ p ) + (_i)"(W)) J] e(o, jf«fl-W 

a-1 M 

xn e (^ a ) H(j ' a;P)+ ^' 6;P) IT e(a,j) P(aj;P)+H(6j;P) II e(j,£;)^ fc;P ), (D.10) 

j=l j=6+l l<j<k<M;j,k^a,b 

and we obtain 

Ajif(P) + A M ((a6)P) = for any P E S M - (D.ll) 

Here we note the following: H(b,a;P) = when H(a,b;P) = 1; H(b,a;P) = 1 when 
H(a,b;P) = 0. 
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Following the discussion in the appendix, we can show the Pauli principle 
of the Bethe ansatz also for the XXZ model; we redefine e(j, k) by e(j, k) = 
sinh(fj — v k + 2r/)/sinh(fj — v t — 2rj), where 77 is related to A in eq. ( p.l|) by A = cosh 277. 

E. APPENDIX: RIGOROUS DERIVATION OF THE COORDINATES BETHE 

ANSATZ 

A. Secular equations of the XXZ model 

The coordinate Bethe ansatz was introduced by Bethe for the one- dimensional XXX 
model in Ref . Q . In this appendix, we derive rigorously some sets of sufficient conditions for 
a vector to be an eigenvector of of the XXZ model under the twisted boundary conditions. 
The derivation should be useful for discussing singular eigenvectors of the model such as 
shown in Refs. [J31-33] . We note that when the twisting parameter is zero: $ = 0, the 



twisted boundary conditions reduces into the periodic boundary conditions, and also that 
when A = 1 the XXZ Hamiltonian (15. If) becomes the XXX Hamiltonian. 



Let us consider the action of the XXZ Hamiltonian Tixxz (5J-) on any given vector. 
We recall that the symbol \M) in eq. ( |2.3| ) denotes a vector with M down-spins where the 
amplitude g(xi,x 2 , ■ ■ ■ , xm) is given by any function. Then, the action of Tixxz on the 
vector \M) can be calculated rigorously. The result is given in the following 

l \ i r M J^ 

Hxxz- JA(M--) \M) = -- Y, EE g{x u ...,xj + 8,...,x M ) 

4 ' Z l<xi<-<x M <L lj=ls=±l 

M-\ . iJ+i j 'i +1 , 

- Y ^+i,x j+ i(g(xi, . . . ,x]^x], . . . ,x M ) +g{x 1 ,...,x j + l,Xj + l,...,x M ) -2Ag(xi,...,x M )) 

-8x 1 ,i8 XMtL (g(0,x 2 , ■ ■ ■ ,x M ) + g(xi,...,x M -i,L + 1) - 2Ag(l,x 2 , ■ ■ ■ ,x M -i,L)) > x J| cr Xh \0) 

J fe=i 

-- J2 \ {-9(0, x 1} ..., Xm-i) + g(xi, . . • , Xm-i, L) e" J * ) af JJ a~ k |0) 

Kx 1 <---<x M _ 1 <L K k=l 



+ (g(l, xi,..., xm-i) - g{xi, . . . , xm-i, L + 1) e" J * ) o~i n a~ h |0) \ 

fc=i J 



(E.l) 
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Here, we have assumed the twisted boundary conditions for the spin operators: & L+1 = 
e ±l ®af] cr£ +1 = of, while any boundary conditions have been assigned on the function 

g(xi,...,x M ). 

Let us discuss sufficient conditions for vector \M) to be an eigenvector of the XXZ 
Hamiltonian, explicitly. Considering the last part of eq. ( |E.1| ) we have the twisted boundary 
conditions on the function g(xi, ■ • • , %): 

g(Q, xi, . . . , x M -i) = g{xi, • • • , x M -u L)e~ l ® for 1 < x x < ■ ■ - < x M -\ < L , 

g(l, xi, . . . , x M -i) = g{xi, ..., x M -i, L + l)e~ J<J> for 1 < x x < - - ■ < x M -i < L . (E.2) 

Considering the second part of eq. ([E.l|) , we have the following conditions. 

gyxi, . . . , Xj, Xj, . . . , xm) + g{xi, ■ ■ ■ , Xj + 1, Xj + 1, . . . , xm) 
-2Ag(x!, ...,xj,Xj + 1,..., x M ) = 
for 1 < xi < ■ ■ ■ < Xj, Xj + 1 < Xj + 2 < ■ ■ • < xm < L and Xj+i — Xj + 1 , (E-3) 

where j is given by j — 1, . . . , M — 1. And we also have 

g(0, x 2 , . . ■ , x M -i, L) + g(l, x 2 , . . . , x M -i, L+l)- 2Ag(l, x 2 , . . . , x M -i, L) = 

for 1 < x 2 < ■ ■ ■ < x M -i < L . (E.4) 



Under the twisted boundary conditions ( |E.2|) , the conditions ( |E.4|) correspond to the special 
cases of the conditions (|E.3|) , where j = M — 1 and xm-i — L, or j — 1 and xi = 0. 

Let us now assume that the function g(xi, . . . , xm) is given by f(x\, . . . , %) defined by 
a general linear combination of the planewave-type solutions 

M 

f(x t , . . . , x M ) = Yl B ( P ) ex P(^ J2 kpjXj) . (E.5) 

Pes M i=i 

Here fci, . . . , Um are free parameters, Sm denotes the symmetric group on M letters, and 
the amplitudes -B(P)'s are arbitrary. The amplitudes B(P) y s in ( |E.5| ) are M! independent 
parameters, and we shall determine them so that the function f(x\, . . . ,%) satisfies the 
conditions ( |E.2|) , ( |E.3|) and ( |E.4|) . The function f(xi, . . . ,%) has the following property 
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M 

/ j \1\X\, • • • ) Xj—\, Xj ~ 1; X j + l-i ■ ■ ■ i X M) + J {X±, • • • 1 X j-ll X j + 1) ^i+1) • ■ ■ J X M)) 

3=1 

' M \ 

y^ 2 cos kj f(xi, . . . , %) , for 1 < x\ < ■ ■ ■ < xm < L . (E.6) 

\i=i / 

In fact, we see explicitly 

M 

Y. (/( x i' • • • , ^- - 1, • • • , xm) + /(an, • • • , Xj + 1, . . . , x M )) 

3=1 

M / A/ A/ 

= JZ 5Z ( 5(P)exp(i^A;p^Xf - ik Pj ) + B(P) exp(iJ2k Pi xi + ik Pj ] 
j=ip<z Sm \ ^=1 £ = 1 

_ M M 

= J2 J2 B ( P ) ex P(^ H k Pi x e ) (exp(-ik Pj ) + exp(ik Pj )) 
PeS M j=i l=i 

M \ 

^ 2 cos % \ f(x!,...,x M )- (E.7) 

Let us summarize the discussion given in the above. Assuming that the function 
g(x\, . . . ,Xm) is given by f(xi, . . . ,Xm) defined by eq. ( E.5|) , the vector \M) is an eigen- 
function of the Hamiltonian if the conditions ([E.2|) , ( |E.3| ) and ( |E.4| ) are satisfied. 

We now show that the conditions ( |E.3| ) are satisfied if the following relations hold for 
the amplitudes _B(P)'s: 

B(Q) (l + exp(ik Qj + ik Q{j +i)) - 2A exp(ik QU+1) )^j 
+ B(Qn j ) [I + exp(ik Qj + ik Q ( j+1) ) - 2Aexp(ik Qj )j = 

for Q G S M and j = 1, . . . , M - 1. (E.8) 

Here the symbol ttj denotes the permutation of j and j + 1: ttj = (j,j + l) for j = 1, . . . , M— 1. 
Explicitly we have 

/ (xi, . . . , Xj—i, Xj, Xj, Xj+2-, • • • , Xm) + J\Xii • • • , Xj—ii x j + 1, Xj + 1, Xj+2, • • • , Xm) 
— 2A/(xi, . . . , Xj—i, Xj, Xj + 1, Xj+2, • • • , xm) 
= ^ I B(P) exp J3 i/ep^x* + ikpjXj + ik P{j+l) Xj 

P&Sm L \£=l;ift,j+l / 

(M 
y^ ikpiXi + ikpjixj + 1) + ikp(j + i)(xj + 1) 
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—2AB(P) exp J^ ikpeX£ + ikpjXj + ikp(j + i)(xj + 1) \ > 

\e=i;e^j,j+i ' / 

/ M 

Y^ B(P)exp\ Yl ikpexe + i(k Pj + k P{j+1) )xj 

Pes M \t=l; tyj,j+l 



(l + exp(ikpj + ikp(j+i)) - 2Aexp(ikp(j+i)) 

1 / M 

9 E ex P E ik Qe x t + K k Qj + k QU+i)) 

z qgs m \e=i- e&,j+i 

x\b(Q) (l + exp(ik Qj + ik Q{j+1) ) - 2Aexp(ik Q{j+1)j 

+B(Q7Tj) (1 + exp(ik Qj + ik Q ( j+1) ) - 2Aexp(ik Qj )) \. (E.9) 



Here, we have made use of the following relation 

E np) = \ E f{Q) + \ E no*i). (E.10) 

P&S M Z Q&Sm Q&Sm 

There are (M — I) Ml/ 2 independent relations in eq. ( |E.8| ). If they hold, then the terms 
involving B(Q)'s and S(Q7Tj)'s in RHS of eq. flE.9Q vanish. Thus, LHS of (|E.9| ) becomes 



zero, and the conditions ( E.3|) are satisfied for the case of j. Thus, we have shown that the 
relations ( E.8|) are sufficient for the conditions ( E.3|) . Hereafter, we shall call the conditions 
( E.8|) the vanishing conditions. 

Let us discuss the number W of independent relations given by eqs. ( E.3|) and (|E.4j ). It 
is given by the number of configurations where x\, . . . , xm satisfy the conditions: 1 < X\ < 
■ ■ ■ < Xj, Xj + 1 < Xj + 2 < ■ ■ ■ < xm < L and Xj+i — Xj + 1 for j — 1, . . . , M — 1 and also 



those of ( |E.4| ). The number W is given by 



We recall that the number V of the relations of ( E.8|) is given by Ml. 

Let us now consider the ratio W/V. We recall that the vanishing conditions ( E.8| ) are 
sufficient for the conditions ( E.3| ) and ( |E.4| ) to hold. If the ratio W/V is larger than 1, then 
the vanishing conditions ( |L.8| ) are also necessary conditions for ( |E.3| ) and ( |E.4| ). Here we 



recall that the variables hi, . . . , ku are assumed to be free parameters in this subsection. 
Let us calculate the ratio W/V, explicitly. It is given by 
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W/V= 2L{L ~ 2)] (E.12) 

1 (L - M)\M\(M - 1)\ y ' 

For example, we have W/V = L for M = 2; W/V = L{L - 2)/6 for M = 3; W/V = 

L(L — 2)(L — 3)/72 for M = 4. Thus, when L is large enough with respect to the number 

of down-spins M, then the ratio W/V is larger than one. When W/V > 1, the vanishing 

conditions (|E.8|) are necessary for the relations ( |E.3|) and ( |E.4|) . However, we should note 



that the ratio W/V is not always larger than 1. For instance, when L = 16 and M = 8, 
we have W/V = 143/420. For the half-filling case, we have W/V < 1 for M = L/2 > 8, 
in general. If W/V < 1, then the vanishing conditions ( |E.8| ) for the amplitudes are not 
necessary for the relations ( |E.3| ) and ( |E.4|) . 

Let us make a conclusion of Appendix E.A. Assuming that the function g(x\, . . . , %) is 
given by f(xx, . . . ,xm) defined by eq. ( |E.5| ), the vector \M) is an eigenvector of the XXZ 
Hamiltonian, if the vanishing conditions ( [F.8| ) and the twisted boundary conditions ( |F.2| ) 
hold. 

We give remarks. Let us consider the case of |A| < 1, where A = cosh 27/. When all the 

momenta are generic, or kgj ^ ±2|?7| (mod27r) or A;q(j+i) ^ rJ=2 {77 1 (mod27r), the relations 

( |E.8|) can be expressed as follows 

B{Q) = 1 + exp(ik Qj + ik Q{j+l) ) - 2A exp(ik Qj ) 

B(Q-Kj) 1 + exp(ik Qj + ik Q{j+l) ) - 2A exp(iA;Q (i+ i)) ' 

for QeS M and j = 1, . . . , M - 1. (E.13) 



For the singular solutions of the XXZ model such as discussed in Ref. H3T|-p3|], some of them 
satisfy the sufficient conditions ( |E.8|) and ( |E.25| ), but their amplitudes S(P)'s do not satisfy 
the factorization property ( |E.13j ). Some details will be discussed elsewhere. 

B. Amplitudes of the Bethe ansatz wavefunction 



We show that the amplitudes Am(PYs defined by eq. ( |1.3| ) satisfy the vanishing con- 
ditions ( |E.8j ), where the parameters ki, . . . , ku are generic. Let us introduce the following 
notation 

46 



F jt (P) = exp[i(k Pj + kpt)] + 1 - 2A exp(ik Pj ) (E.14) 

Then, the expression (|1.3|) of the amplitude Am{P) is given by 



A M {p)=e{p) n §4t' for peSm - ( e - 15 ) 

i<j<e<M ^3^) 



Here e denotes the unit element of the permutation group Sm, and we put C — 1 in (|L 

Let us calculate the ratio of A M (Qir a ) and Am((3), explicitly. Here a is taken to be an 
integer satisfying 1 < a < M — 1. We note that integers j and £ satisfy the condition: 
l<j<£<Mva ( |K.15| ). We consider the four cases: (1) j = a and £ = a + 1; (2) j — a 

and £ > a + 1; (3) j = a + 1 and £ > a + 1; (4) j, £ 7^ a, a + 1. Then we have 

A M (Q)Y[F je (e) = e(Q)F aa+1 (Q) J[ F a£ (Q) J[ F a+U (Q) J[ F j£ (Q) (E.16) 

j<£ £>a+l £>a+l l<j<£<M;j,e=£a,a+l 

When P = QiTa, we have 

A M {Q^a)\{F 3 ^)=<Q^a)Faa + l{Q^a) \{ F ai {Q<K a ) \[ F a +u(Q^a) 

j<£ £>a+l £>a+l 

x J] F jt (Qir a ) (E.17) 

j<£;j,ey£a,a+l 

Through an explicit calculation, we have 

F aa+ i(Qn- a ) = F a+lta (Q) , F a+la (Qn a ) = F a , a+1 (Q) 

FaiiQ^a) = Fa+iAQ) , F a+u {Qix a ) = F a>i (Q) (E.18) 



Thus, the ratio is given by 



A M (Q) _ , ^ F ata+1 (Q) 



A M (Q^a) F a+lA (Q) 

It is nothing but the relation ([E.13|) , which is equivalent to the vanishing conditions ( E.8| ) 
for the case of a. Here we recall that all the momenta are given generic in Appendix E.B. 

Let us now discuss the uniqueness or well-definedness of the amplitudes v4 M (P)'s. First 
we note that any permutation P can be written in terms of a product of generators 7Tj's. For 
example, cyclic permutation (213) corresponds to (23) (12) = 7r 2 7ri. Thus, we can calculate 
Am{P) by using the relations ( E.13|) . For example, let us take P = (213). We have 
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A(7r 2 7Ti) A(7T 2 ) ,,, 

A(7F27ri) = A(vr 2 ) A(e) A(6) 
= ^21(^2) ^32 (e) 
i^fo) F 23 (e) l ' } 

However, different products of generators can correspond to the same permutation. For 
instance, we have (213) = (23)(12) = (12)(23)(12)(23). The amplitude A(ir 2 Tn) should be 
equivalent to v4(7ri7r 2 7Ti7r 2 ). 

In fact, we can prove the uniqueness of amplitudes Am(P)'s, explicitly. We note that 
the defining relations |37[ of the symmetric group Sm given by the following: 



tt| = 1 for j = l,...,M-l, (E.21) 

-Kj-Kj+iTij = irj+iTtjTtj+i for j — 1, . . . , M — 2 (E.22) 

Thus, for a given permutation P, any given two products of generators expressing the same 
P can be transformed into one another, by using the defining relations given in eqs. ( E.21|) 



and ( E.22|) . Therefore, the uniqueness of the amplitudes is proven if we show that the 
amplitudes Am(-P)'s satisfy the following relations: 



A M (Q) A M (Qnj] 



Am(Qkj) A M ({Qirj)irj) 

A M {Q) A M (Qnj) A M (Qir j ir j+1 



1 for j = 1, . . . , M - 1 , (E.23) 



A M {Q-Kj) A M ((Q7rj)7r j+1 ) A M (Q7r j 7r j+1 7T j ) 

. for j = 1, . . . , M - 2 . (E.24) 



A M (Q) A M (QiTj +1 ) A M (Qn j+ i7Vj) 



Am(Qkj+i) A M ((Q^j+i)^j) Am{{Qt^j+x'Kj) / Kj+x) 
Here Q G Sm- in fact, it is easy to check the relations (E.23|) and ( E.24| ). 



C. Derivation of the Bethe ansatz equations 

The twisted boundary conditions ( |E.2[ ) for the wavefunction are given by the following 
f(xi, . . . , x M ) = exp(-i$)/(x 2 , . . . , x M , Xi + L) for < x x < ■ ■ ■ < x M < L (E.25) 
In terms of the amplitudes 5(P)'s, RHS of (E.25|) is given by 
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J2 B(P) exp (-i$ + i(k P1 x 2 H k P ( M -i)X M + k PM (xi + N)) 



Pes M 



J2 B(P) exp (-2$ + ik PM N + i(k PM Xi + k P1 x 2 H k P(M _ x) x M \ 



PeS A _ 

M 



= ^ B(Q(12---M))e^p\-i^ + ik Q1 N + iY l k Qj x j \. (E.26) 

Qg5 m V i=i / 

Here we note that for Q(12 • • • M) = P we have 

Q1 = PM, Q2 = P1,..., QM = P(M-l). (E.27) 



In order for RHS (|E.26| ) to be equivalent to LHS of (|E.25| ), we have 



exp(ife gi iV) = exp(i^) B *®\ M)) , for QeS M . (E.28) 



Let us now derive the Bethe ansatz equations from the eqs. ( |E.28| ). Here, we assume 



that the amplitudes P(P)'s are given by Am{P)'s satisfying the vanishing conditions ( |E.8| ) 
(or ( |E.13| )). Then , we have the following 

A(Q) A(Q) Ajm) A M (Q7T 1 ---7r M -2) 

A(Q(12 ■ ■ ■ M)) A(Q<kx) ' AiQir^) ' ' ' A m ((tt 1 ■ ■ • tt m _ 2 ) 7r M -i) 

= f-i)^- 1 T\ ( l + ex P( i ( fc Qi + k Qm) - 2Aexp(ik Q1 ) \ 
rn=2 V 1 + exp(i(/cQi + k Qm )) - 2A exp(ik Qm ) J 

(E.29) 

Here we note that Qtti ■ ■ ■ 7iy_ir = Ql and Qir\ ■ ■ ■ 7r r _i(r + 1) = Q(r + 1). Thus, we have 

cti \ t -\\M-i ok\ TT ( l + exp(i(/c Q i + fc Qm ) -2Aexp(i/c i) \ 

exp(«Lfc Q i) = (-1) exp(?$) {{ — — - — — - — - — — — — - 

m=2 V 1 + exp(«(fc Q1 + fcg m )) - 2Aexp(ik Qm )J 

for Q G S M - (E.30) 

Let us write Ql by j. Then, Q2, . . . , QM are given by all the integers from 1 to M except 
j. We may write Qm by £ which runs from 1 to M except j. Thus, we obtain the standard 
form of the Bethe ansatz equations. 

Finally, we summarize the rigorous formulation of the coordinate Bethe ansatz. Assuming 
that the function g(x\, . . . ,Xm) is given by f(xi, . . . ,Xm) defined by eq. ( |E.5| ), we have 
shown in Appendix E.A that the vector \M) is an eigenvector of the XXZ Hamiltonian, if 
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the vanishing conditions (|E.8|) and the twisted boundary conditions ( |E.2|) hold. In Appendix 
E.B, we have shown that the amplitudes Ajy(P)'s defined by eq. ( |1.3| ) are well-defined and 
also that they indeed satisfy the vanishing conditions ( |E.8| ). In Appendix E.C, we have shown 
that the twisted boundary conditions ( |E.2|) are satisfied when the Bethe ansatz equations 
( |E.3U| ) hold. Therefore, if k%, . . . , ku satisfy the Bethe ansatz equations, with the amplitudes 



Am"(-P)'s constructed by the momenta via eq. (|1.3| ), the vector \M) becomes an eigenvector 
of the XXZ model under the twisted boundary conditions. 
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